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Abstract 

We construct fractional Brownian motion (fBm), sub- fractional 
Brownian motion (sub- fBm) , negative sub- fractional Brownian motion 
(nsfBm) and the odd part of fBm in the sense of Dzhaparidze and 
van Zanten (2004) by means of limiting procedures applied to some 
particle systems. These processes are obtained for full ranges of Hurst 
parameter. Particle picture interpretations of sub- fBm and nsfBm were 
known earlier (using a different approach) for narrow ranges of param- 
eters; the odd part of fBm process had not been given any physical 
interpretation at all. 

Our approach consists in representing these processes as (AT(1), l[o,ti), 

l[ ,t] - l[-t,o])) (-X"(l)»l[-t,t])> respectively, where X(l) is an 
(extended) <S'-random variable obtained as the fluctuation limit of ei- 
ther empirical process or the occupation time process of an appropriate 
particle system. In fact, our construction is more general, permitting 
to obtain some new Gaussian processes, as well as multidimensional 
random fields. In particular, we generalize and presumably simplify 
some results by Hambly and Jones (2007). We also obtain a new class 
of S'- valued density processes, containing as a particular case the den- 
sity process of Martin-L6f (1976). 
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1 Introduction 



The main objective of the present paper is to give a particle picture in- 
terpretation of the continuous centered Gaussian processes £ H = (£t*)t>o, 
Q H = (C/0t>o an d V H = ( r ??)t>o with covariance functions, for < s < t, 

Egg =\{s 2H + t 2H - \t - s\ 2H ), H e (0, 1) (1.1) 

E(? Cf =(1 - H) (s 2H + t 2H - + tf H + |t - s\ 2H 

He (0,2) (1.2) 
r/f =( S + t) 2 " - |t - s| 2 " , € (0, 1). (1.3) 

£ H is the well-known fractional Brownian motion (fBm), £ , for i7 < 1, is 
the sub- fractional Brownian motion (subfBm, see e.g. [3]), and for H > 1 
is the so-called negative sub-fractional Brownian motion (nsfBm [7]), r] H is 
a counterpart of the subfBm in the sense that both processes occur in the 
decomposition of fBm studied in [13J. It is called the odd part of fBm in 
that paper. It is also related to nsfBm, namely, = K r]?~ 1 ds for a 
constant K (see [7]). It is known that the closures of the intervals written 
above are maximal ranges of the Hurst parameters H. In each case, the 
extreme points of these intervals correspond to trivial processes, which will 
be excluded from our considerations. 

There are many models, related to particle systems or not, which lead to 
fBm (see, e.g., @],[26], p2], [H]), especially, for H > ~. Sub- fBm and nsfBm 
appeared in a natural way in connection with occupation time fluctuations 
of particle systems ([12], [2], [3], [8], [9], we give some more details below) 
and this was the reason for us to study them in [3] and [Tj. Sub- fractional 
Brownian motion has gained independent interest and was investigated also 
by other authors (e.g. [JJ, [27], [ID], [28]). However, the particle models 
studied earlier led to these processes for narrow ranges of parameters only; 
namely, H G [|, 1) for sub-fBm and H G (1, |] for nsfBm. It seemed natural, 
and interesting, to ask if there exist other particle models that would permit 
to obtain these processes for the whole ranges of parameters. For example, 
in the case of sub-fBm, the difficulty of interpretation seemed to follow 
from the fact that for H < ^ increments of this process on non-overlapping 
intervals are negatively correlated. As far as we know, the process rj H has 
not been given any physical interpretation at all. 

In this paper we use "the white-noise approach", i.e., our starting point is 
the well-known fact that the standard Brownian motion can be represented 
as (X, l[o,t]) 5 where X is the white noise (see, e.g. [2TJ). A similar con- 
struction permits to obtain processes defined in (jl.ip — (|1.3p for all possible 
H. It also yields new interpretations of sub-fBm and nsfBm for parameters 
obtained before. 

The second objective, related to the first one, is to introduce a new 
class of density processes in 5', and associate it to a class of real Gaussian 
processes, containing as a special case the process studied by Hambly and 
Jones (P5], HU). 
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All the models discussed in this paper are based on the following particle 
system in R: at time t = positions of the particles are determined by a 
point measure v (in general random), which is in some sense homogeneous; 
then they are evolving independently according to the standard, symmetric 
a-stable Levy motion (a £ (0, 2]). In some models they additionally undergo 
a critical binary branching. The evolution of the system is described by the 
empirical process N = (Nt)t>o, where ^(^4) is the number of particles in 
the set A C M at time t. In general, the system depends on a parameter 
T which will tend to infinity. The corresponding empirical process will be 
denoted by N T . 

We consider three classes of models. The first one gives a particle picture 
interpretation of £ , ( H and r] H for H < | in a somewhat circuitous way 
(two limit passages), but, on the other hand, it is quite general and leads 
to a new class of the so called density processes. Moreover, we obtain a 
convergence result of Hambly and Jones |16j . |17j . as a very particular case. 
The initial configuration is given by a measure vt such that as T — > oo the 
density of particles increases. 

We define 

X T{t) = ?l^L. (1.4) 



If vt is a Poisson random measure with intensity measure TX, where A 
is the Lebesgue measure, it is well known that Xt, regarded as 5'-valued 
processes (5' is the space of tempered distributions) converge in law to 
the so-called density process, which is a continuous centered Gaussian S'- 
valued process of Ornstein-Uhlenbeck type (see e.g. [23] , [IB] , [15] ) • Our 
first result is a generalization of this fact (in the non-branching case). We 
consider a class of initial configurations vj- which includes as particular cases 
both homogeneous Poisson random measure and the deterministic measure 
Y2jezfi±- We show convergence of finite dimensional distributions in S' to 
a certain new density process X (Theorem I2.ip . Moreover, we show that 
this convergence holds for a class of test functions wider than S (the space 
of smooth rapidly decreasing functions). We then study real processes of 
the form 

((X(l),1> t ))t>a, (1.5) 

where {X(l),ip) is defined by a natural extension. 
The most important examples are 

^ =l[o,t] (1-6) 
ik =l[ ,t] " l[-t,o] (1-7) 
i>t =l[_t,t] (1-8) 

Note that the space parameter plays now the role of "time" . For ipt of the 
form (|1.6p we obtain a generalization of a convergence result by Hambly and 
Jones ([E], [E]) who took a diffferent approach and did not use explicitely 
the high density of the system (see Remark l2,4l (b)). Note that the processes 
(Ni,ipt) m each case (jl,6p — (|1.8p have a clear physical interpretation. For 
(jl.7p the interpretation of (jl.4|) is particularly nice, since E (Nf \ipt) = if 
vt is symmetric. 
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Next, assuming that the initial configuration is deterministic or nearly 
deterministic, we make a simple second passage to the limit. We speed up the 
"time" in (|1.5p (in fact, shrinking the space), i.e. we consider (X(l), tpTtji or 
more generally, (X (1) , 92(7^)) and normalize appropriately. In the limit, as 
T — > 00, for ipt of the form (jTTHJ) we get easily a fractional Brownian motion 
(this result was proved by Hambly and Jones [16] ) , for (jl.7p a sub-fractional 
Brownian motion, and for (|1.8p the process rj H defined in (|1.3p . with H < ^ 
in each case (Propositions 12.51 12 .6[) . In this way we obtain the desired 
particle picture interpretation of the sub-fBm for small Hurst parameters. 
We remark that this procedure for truly random initial measures does not 
lead to interesting results. 

It is worthwhile to note that writing the processes in the form (|1.5p 
permits to derive easily their properties, such as stationarity of increments, 
negative correlatedeness of increments, long range dependence. Some of 
these properties are preserved after the second passage to the limit. 

In the second model, giving the particle picture interpretation of £ H , £ , 
r] H for H < ^ more directly, we consider an initial configuration determined 
by a fixed measure v, which is deterministic, or almost deterministic. We 
assume a < 1; again, there is no branching. We carry out the same space 
transformation as in the second passage to the limit in the previous model 
and define 

{ZT ^AjhM^0iM^l, (1 . 9) 

We prove convergence in law in S' to a centered Gaussian ^'-variable. Again, 
this result can be extended to test functions ipt of the forms (jl.6p - (|1.8p . yield- 
ing convergence of {Zt,iP-} in the sense of finite dimensional distributions 
to fBm, sub-fBm, r) H , respectively, with H < ^ (Theorem 12. 9p . 

The third class of models is related to occupation time fluctuations. We 
consider a particle system with homogeneous Poisson initial condition with 
or without high density, with or without branching and define 

Y T (t) = jr ( jT* Njds - J™ EN^ds^j . (1.10) 

where Ft is a deterministic norming. 

In [5] and [8] , for "large" a (a > 1 in the case without branching, and a > 
^ in the case with branching), we showed the weak functional convergence 
of Yt in S' , as T — > 00, to a limit process Y which is of the form K-dX, where 
i9 is a centered Gaussian process. Depending on the specific parameters of 
the model, i9 is 

• a fractional Brownian motion with H 6 (^, |], 

• a sub fractional Brownian motion with H G (^, 1), 

• a negative subfr actional Brownian motion with H £ (1, j). 

Considering a branching system in equilibrium Milos in [23] obtained a frac- 
tional Brownian motion for H G (^, 1). On the other hand, for "small" a 
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the limit Y also exists ([5]) but is of a completely different form. It is truly 
5'-valued process with independent increments. 

In this paper, in the non-branching case with a < 1, we consider (Yt(1), ipt) 
for ipt of the form (jl.6p - (|1.8p . similarly as for density processes, and in the 
limit we obtain fBm, sub-fBm and the process rf H with H G 1). The fact 
that the formal expression (Y(l), l[o,t]) leads to a fBm was observed in [22j . 
Moreover, for a > 1, considering (YT,ipt) with ipt of the form (|1.7p we get 
a nsfBm with if G (1, |). In this case the norming Ft is different from the 
one in [3], due to the fact that J R ipt(x)dx = 0. 

So we see that as a changes from "small" to "large" values, there occurrs 
a sort of phase transition in the sense that a complicated time structure 
and simple space structure changes into a simple temporal structure and 
complicated spatial structure. Moreover, this "complicated" structure in 
both cases is related to fractional Brownian motion. 

It remains to find a particle picture yielding a negative sub-fractional 
Brownian motion for the missing range of parameters H G [§ 5 2). This is 
achieved by considering a high density branching particle system with high 
branching intensity (Theorem l2.15p . 

We remark that the argument of this paper can be carried out in the 
multidimensional case, for example we can consider ipti,...,t d = l[o,ti]x...x[o,t (i b 
obtaining some Gaussian random fields. 

It seems more difficult to try to extend the approach of this paper to 
the systems with infinite variance branching. It is worthwhile to notice that 
taking formally (X(l), l[o,t]), where X is the stable 5'-valued limit process 
in Theorem 2.1 in [6], one obtains a linear fractional stable motion with 
parameters a = 1, b = —1, ^ < H < 1 defined in |25| (Definition 7.4.1). One 
would obtain an interesting and probably new interpretation of this process 
if this procedure were justified by a suitable limit theorem. 

The following notation is used in the paper: 
A: Lebesgue measure; 

S: space of C°° rapidly decreasing functions on R; 
S': space of tempered distributions (topological dual of S); 
(•, •): duality, in particular on S' x S and f) = f fdfj,; 
=^-: weak convergence of finite dimensional distributions of processes in ap- 
propriate space; 

Pt(x): transition probability density of the standard symmetric a-stable 
Levy process in R; 

Tt- semigroup determined by pt, i.e. Tt^p = Pt* </>; 
= Jr e ixy ^{y)dy. 

Generic constants are denoted by C, Cj with possible dependencies in paren- 
thesis or in subscripts. 

In Section 2 we describe the particle systems, formulate the results and 
discuss them. Section 3 contains proofs. 
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2 Results 



2.1 Density processes 

We consider the following particle system introduced in [9]. 

Let 9 be a non-negative integer- valued random variable with distribution 

P(9 = k) =p k , k = 0,1,2, ... , (2.1) 

and such that E6 3 < oo. Let 6j,j G Z, be independent copies of 0, and 
for each j € Z and = 1,2,..., let p^. = (fP k v . . . , f? k k ) be a random 

vector with values in [j,j + l) k . We assume that (6j, (pL)k=i,2, ■ ■ •)> 3 £ Z, 
are independent. Given T > 0, these objects determine a random point 
measure vt on R in the following way: For each j, 6j is the number of points 
in the interval [^,^f-) ; and for each k, if #j = k, the positions of those 
points are determined by ^pj,- In other words, 

j £ Z n=l 

where ^ 

Kj,n,T = yPj^n, (2-3) 

and 5 a is the Dirac measure at a E R. 

Observe that both the deterministic measure 

and the Poisson random measure with intensity TX are of the form (|2,2p . 

Fix a G (0,2] and assume that at the initial time t = there is a 
collection of particles in R with positions determined by a measure z/t of the 
form (|2.2p . As time evolves, these particles move independently according 
to the symmetric a-stable Levy process. 

We define a signed measure- valued process Xt by (jl.4p . for which we 
have the following result. 

Theorem 2.1. =4> X in S' as T — > oo, where X is a centered Gaussian 
S' -valued process with covariance functional 

E(X(s),(p) (X(t),<ip)=E6 [ ip(x)T\t- s \^(x)dx+(Y&r6-E6) [ T s ip(x)T t ^(x)dx. 

(2.5) 

Remark 2.2. (a) The limit process X will be called density process. Note 
that if vt is a Poisson random measure with intensity TX or, more generally, 
if 9 is such that Ed = Var#, then X is the classical density process (e.g. 
E3]). 
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(b) The same result is true if the a-stable motion is replaced by a more 
general process with the covariance functional of the limit written as 

E9 [ Eip{x + $ s )if){x + '& t )dx + (Var0-£0) f E<p(x + ti s )Eif)(x + ti t )dx. 
Jr Jr. 

(2.6) 

The only condition on ■d is that there exists m > such that (|2.6j) is finite 
if ip and if) are replaced by (f> m (x) = 1+ y*n ■ 

(c) There is no problem to extend this result to M. d . In the definition of vt the 
interval [j, j + 1) should be replaced by the cube \j\ , j± + 1) X . . . X [jd, jd + 1) 
and in (|1.4p the normalization is v / T^. The limit is an ^(R^-valued process 
with covariance of the form (|2.6p . 

We want to extend (X, ^) for some ijj's which are not necessarily in S. 
Let \l/ be the class of bounded piecewise continuous functions with compact 
support. From (|2.5p it is clear that for ip £ \1/ we can define (X, ^) by L 2 - 
approximation. Let \& = (4>t)t>o be a family of functions from Typically, 

are of the form ()1.6p - ()1.8p . We will study (centered Gaussian) processes 
of the form 

gf = (X(l),i/) t ) t > . (2.7) 
By (|2.5p . the covariance function of is 

Egfgf =E6 f ip s (x)if) t (x)dx + (Var0 - £0) / i) s (x)T 2 Mx)dx. (2.8) 

The next proposition shows that this extension is compatible with The- 
orem 

Proposition 2.3. 

((X T (l),if) t )) t > =» asT^oo. 

Remark 2.4. (a) This proposition gives a particle picture interpretation of 
the process g^ . 

(b) For deterministic 9 and of the form (|1.6p . the process is, up to 
a constant, the same as the process G defined in Proposition 4.2 of |16j . 
(see also [T7j) for parameter c = 1. Moreover, the convergence result of that 
proposition follows directly from Proposition [2]3] by a standard conditioning. 
(The extension to general c is also immediate.) We stress that the proof in 
|16j is completely different, not exhibiting the high density of the system. 

The covariance (|2.8p may look complicated but it does permit to deduce 
some properties of the underlying process (see Proposition 12. 7p . The sit- 
uation simplifies considerably after the second passage to the limit which 
consists in shrinking the space. 

For a function / : M H> 1 and T > we denote 

S T f(x) = /(£). (2.9) 
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We define an 5'-valued random variable Zt by 

Zt,<p) = ^-(XQ),St<p), (2.10) 

where X is given by Theorem 12.11 and Ft is a suitable norming. We have 
the following simple fact. 

Proposition 2.5. (a) If Var 9 = and Ft = \/T^~ a then Zt converges in 
law, as T — ^ oo, to a centered Gaussian S' -valued random variable Z with 
covariance 

E (Z, (p) (Z, ip) = — / (p(x)ip(x)\x\ a dx. (2.11) 

7T ./id 



(b) IfY&iO > and Ft = vT then Zt converges in law in S' to (Var 9) W , 
where W is the white noise. 

From the proof it is clear that in this proposition the functions ip £ S 
can be replaced by ifit £ ^> giving the convergence of finite dimensional dis- 
tributions of the processes -p-Q. T * := -p- {X(l), Srifi-} (cf. (|2.7|) ) provided 

2 

that f R ifit(x) \x\ a dx < oo in case (a). In particular, if $ has one of the 
forms (fO jl -tfTS ]) . then 

Qt = Qrt I 2 - 12 ) 
and we have the following Proposition: 

Proposition 2.6. Assume Var# = and a < 1. Let \& /iaue one o/ f/ie 

forms (|1.6|) - (|1.8|) . Then the processes ( ^ T \_ a Qrt )* >0 conver 9 e *n f/ie sense 

of finite dimensional distributions to: 

(a) a fractional Brownian motion K^ H if ifit = l[o,tp 

f&j a sub-fractional Brownian motion KC, H if ifit = l[o,t] ~~ ^[-t,ol; 

fe,) £/ie process Kn H if ifit = l[-t,t]> 

with H = in each case. 

We have assumed that Var # = 0, i.e., the initial configuration of particles 
is close to deterministic, since for random 9, in all cases (|1.6p - (jl.8p one 
obtains the Brownian motion. This proposition gives the desired particle 
picture interpretation for sub-fBm with H < ^ as well as an interpretation 
of the process 77 . Part (a) was observed in [16]. 

We close this subsection with a brief discussion of the properties of the 
process 0* and of the iS'-random variable -^(1). 

Proposition 2.7. Fix \P = (ifit)t>o> ifit £ ^> an d let £>* 6e i/ie process 
defined by (|2.7p . and X(l) 6e given by Theorem \2.1\ 

(a) Homogeneity: (X(l), ip) has the same distribution as (X(l), (p(- — a)) for 
a G R ; (p G S, and, for ^ of the form (jl.6p . has stationary increments. 

(b) Correlation of increments: Assume that ifs<t<u<v then ifit—tfi s and 
ifiv — ifiu have disjoint supports; moreover, the function t 1— > ifit(x) is either 
non- decreasing for x G R or this property holds on R + and the functions tfi t 
are odd. Then 
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(i) if Var 9 < E9 then the increments of g^ are negatively correlated; 

(ii) if Var 9 = E9 then g^ has independent increments; 

(Hi) if Var 9 > E9 then the increments of g are positively correlated. 

(c) Long range dependence: Ifipi,<f2 have compact supports then 



lim r 1+a E {Xil),^) (X(l),^ 2 (- ~ t)) =C{Vax9-E9) I Vl (x)dx I <p 2 (x)dx. 

2.13) 



In particular, for ^ of the form (|1.6p 

lim r 1+a E (gf - g*) {g* +T - g* +T ) = C(V a i9-E9)(t-s)(v-u). (2.14) 

(d) Path continuity: If 

^ tl {x)-^ t2 {x)) 2 dx<C{T)\t l -t 2 f , h,t 2 <T, (2.15) 



JR 



for some 2 > j3 > and any T > 0, then has a continuous version; more 
precisely it is locally Holder continuous with exponent < 

All these properties follow easily from (|2.5p . f|2.8j) . We will give a brief 
explanation in the next section. 

Remark 2.8. (a) A similar calculation as in deriving ()2.14p shows that the 
long range dependence rate of increments of g^ with ^ of the form (jl.8p is 
also r~( 1+a \ On the other hand, this rate for ^ of the form (|1.7j) is t~( 2+Q! ) 
due to the fact that in this case ife are odd. 

(b) Properties (a) and (b) are clearly preserved after the second passage to 
the limit. In particular, we thus obtain a simple proof of negative correlat- 
edness of increments of sub- fractional Brownian motion for H < \. This 
fact was proved in [3] but the argument was rather cumbersome. 

(c) The relationship between the sign of Var# — E9 and positive/negative 
correlatedness of increments seems quite interesting and unexpected. In 
particular, we see that independence of incrementss occurs not only in the 
Poisson case. 

(d) In the special case Var 9 = and $ of the form (jl.6p property (|2.14p and 
stationarity of increments of g^ was obtained in [16] (see Remark 12.4( b)). 

(e) In general, the process g^ is not self-similar. Nevertheless, if Ft is reg- 
ularly varying at infinity, then the =4> limit of ^-^£>f^ , if it exists, is a 

self-similar process (cf. (|2.12p and Proposition 12. 6p . 

(f) This proposition can be reformulated for more general particle motions 
and for the multidimensional case. 



2.2 Particle picture for small Hurst parameters; direct ap- 
proach. 

As announced in the Introduction, in this subsection we show how to obtain 
fBm, sub-fBm, r/ H with H < |, employing just one passage to the limit. 

Consider a particle system described in the previous subsection with 
initial configuration determined by a measure v given by (|2.2p with T = 
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1 and deterministic 9. Recall that this means that the initial number of 
particles in each interval [j,j + 1), j G Z, is fixed, non-random (see (|2.1|) ). 



Theorem 2.9. Assume that a < 1 and define Zt by (JT79J) . 

(a) Zt converges in law, as T — >■ oo, in S' to the random variable Z as in 
Proposition \2.5\f a). 

(b) Let (ipt)t>o have one of the forms (|1.6p - (|1.8p . Then 



{{Z T ^ t )) t > Q ^Kd H , as T^oo, 



where 





if ipt = 


l[0,t]i 




< C H 


if ipt = 


l[o,t] - 1 [-t,o] 1 


(2.16) 




if ik = 


l[-t,t]. 





with H 



l-a 



2 ' 

Remark 2.10. Theorem [231 has an immediate extension to the <i-dimensional 
case. The condition a < 1 is then replaced by a < d. 

2.3 Processes resulting from occupation limits 

So far we have obtained the particle picture for the processes £ H , C H and 
rj H with H < If. In this subsection we show how our "white noise" approach 
permits to obtain them for H > ^. 

We consider the same particle system as before, with or without branch- 
ing, we assume that the initial configuration is homegeneous Poisson. We 
define Yp by (jl.lOp with an appropriate Ft- (If the system does not depend 
on T we write N instead of N T ). 

Proposition 2.11. Assume that there is no branching, the initial measure is 
Poisson with intensity A and a < 1. Let Ft = vT and consider ^ = (ipt)t>0; 
i/) t 6*. Then 

«Y T (l),^)) t > =^ as T^oo, 
where $ is a centered Gaussian process with covariance 

Mrf, = - Mx)Mx)rLdx. (2.17) 
vr ./to \x\ 



Corollary 2.12. If ^ has one of the forms (|1.6p - (|1.8p then the limit # has 
the form i? t = Ktif where $ H is as in (f!TL6|) with H = ±±^. 

Remark 2.13. (a) The covariance (|2.17p can be written in a more explicit 
form 

2i-«r(±=s) /• ^( X )^( X ) 



V^ r (f) Vr 2 |ar — y| 
(b) The process i? which appears in Proposition 12.111 can be obtained by 
putting formally (WQ a \l), ipt\ where the S'- valued Wiener process Wq°^ 
is the weak functional limit of Yt (see Theorem 2.1 in [5]). Proposition 12 . 1 ll 
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shows that this formality is justified. The fact that (W^ , (1), l[o,*] / i s a 
fractional Brownian motion was observed in [22]. 

(c) For "large" a, fBm and sub-fBm with large H appeared in the temporal 
structure of the limit of Yt in models with or without branching and with 
Poisson, equilibrium or deterministic initial conditions (see [3], [23], [9]). 
Also nsfBm with H G (1, | ) was obtained in this way from a branching 
Poisson system with high initial density and a > 1 (see [8]). 

In the next proposition we show that the negative sub-fractional Brow- 
nian motion, and with a wider range of parameter H, appears also in the 
spatial structure of the limit. 

Proposition 2.14. Assume that there is no branching, the initial measure 
is Poisson with intensity A, a > 1 and Ft = yT. For tpt of the form (jl.7p 

((Y T (l)^ t )) t > jK( H as T^oc, 

withH = G (1, §]. 

The covariance of Q H has again the form (|2.17|) , where the corresponding 
integral is finite due to the special form of ipt- In fact, this integral is finite 
for a < 3. The restricted range of H (H < |) follows from our model, which 
requires a < 2. To obtain the particle picture interpretation for nsfBm with 
the full range of parameter H we have to use a more complex model. 

Theorem 2.15. Assume that \ < a < |. Consider binary branching par- 
ticle system with branching rate Vt — > oo. Assume that the initial measure 
is Poisson with intensity HtX such that 

lim T 7 V T H~ l = 0. (2.18) 

T^oo 1 

Let F T = ^TV T H T . Then for ifj t of the form (fTTTjl . 

((Y T (l),iPt)) t >o=>K( H as T^oo, 

with H = G (1,2). 

Remark 2.16. (a) The same model with a < \ gives in the limit a sub-fBm 
with | < H < 1. It seems interesting to observe that ct = \ does not give 
the null process (see (|1.2p ). This is due to the fact (easy to verify) that 

ECS <f = C { H) f a-^WK^ -eosW) ^ (2 19) 

JR \X\ 

for all H G (0, 1) U (1, 2). For H = 1, in the limit we get a non-trivial process 
with covariance (|2.19p . 

(b) We think that interesting results could be obtained for models studied 
in this subsection if branching particle systems were considered with infinite 
variance branching mechanism (see the remark at the end of the Introduc- 
tion). This would require, however, different methods of proof. 
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3 Proofs 

3.1 Proof of Theorem [27T1 

Let i?- 7 '™, j € Z, n = 1,2,..., be independent, standard a- stable Levy 
processes (tfj'" = 0), independent of vt- The process Xt denned by (|1 .4[) 
has the form 

X T [t) = -L J2(N^ - EN?*), (3.1) 

where 7V T, - ? is the empirical process of the system which at time t = starts 
from [j/T, (j + 1)/T), i.e., 

<<^>=£pfo,„,r + ^' n ), (3-2) 

n=l 

see (|2.3p . The processes iV T,J , j G Z, are independent, hence to prove the 
claimed convergence we can use the central limit theorem. To this end, we 
show first the convergence of covariances. 

By independence, (|2.ip . (|3,ip . (|3.2p . we have for s < t, ip,ip E S, 

E(XT^)(Xf^) 

/ oo k k 

=t E X> S (E E + ^ n )^(s, m ,T + tff 

jgZ \A;=0 n=lm=l 
oo k co I \ 

- X>#I>(**«,T + 0j n ) E^^E ^rf + D • (3-3) 



k=0 n=l 1=0 m=l 



Denote h^ n (x) = p£ n — x, where [x] is the largest integer < x. Clearly, 

\h k , n (x)\ < 1. (3.4) 



Splitting the first expression on the right hand side of (|3.3p into the 
sum over the diagonal (n = m) and the rest, by independence and Markov 
property of $, we obtain 



E(Xj,ip)(xT^) = J2pk^2lT(k,n) + J2 Pk £ H T (k, 



n.m) 



k=0 n=l k=0 n,m=l 



+ ^PkPi ^2 III T (k,n;m,l), (3.5) 



k,l=0 n,m=l 

where (substituting x' = x/T) 
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I T (k,n)= f ET s (<p%- s ^)(x+ hk ' n ^ x) )dx, (3.6) 
Jr 1 

f hk n (Tx) hk m (Tx) 
II T {k,n,m) = ET s <f{x-\ '— )T t ip{x H ■ — )dx, (3.7) 

JR 1 1 



ttt (i n I' t?i- I , h k>n(Tx) hi, m (Tx) 
III T (k,n;m,l) = / ET s (p[x-\ '— )ET t ip{x + 

Jr 1 



)dx. (3i 



By (|3.4p . the terms under the integrals in (|3.6[) - (|3.8p converge pointwise 
as T — > oo; on the other hand, if we denote 4>(x) = j^z, then it is easy to 
see that 



T tl (<Pi% 2 (<P2...Tt n <p n )...)(x + y) < CT tl {4>T t2 {<i>...Tt n 4>)..-)(x) (3.9) 
for \y\ < 1, x G R, n = 1, 2, . . . , provided that < Ci</>. Hence 



lim I T (k,n) = / T s ((pT t - s ip)(x)dx = / ip{x)Tt- s tp{x)dx, 
lim IlT(k,n,m) = lim IIlT(k,n;m,l) = / T s <p{x)Tt^{x)dx. 

T— s>oo T— >oo Jjg 

By ([33]), this shows that E (Xj , ip) (X^ , ip) converges to the right hand side 
of (|23D . 

It remains to prove that a i{^T(ti),^>i) converges in law, as T — > oo, 
to a Gaussian variable for any a\, . . . , a m G M, t\, . . . ,t m > 0, (fx, ... , <^ m G 

<S, m = 1, 2, To this end, by (|3.ip . we apply the Lyapunov criterion. It 

is easy to see that the proof will be completed if we show that for any t > 0, 
G 5, > 0, 

lim Arfop) = 0, (3.10) 

where 

ArC*. *>) = E ^fj^ = ^EX> B (E »(w + 

jSZ j& k=0 n=l 

(3.11) 

(see (J321), (EU, J23J). 

Using the obvious inequality [a\ + ■ ■ ■ + a^) 3 < fc 2 (a 3 + • • • + a|) for 
Oi, • • • , a,k > 0, and employing again the functions h^ n we obtain 



Mt,v)<TLjtpk*?jt, I r^\^ + j ^^)dx<^=Ee^ f 4>(x)dx 

(3.12) 

by (J3I1D and (|33g), hence (I3TTTJD follows. □ 
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3.2 Proof of Proposition EO 

Looking at the proof of Theorem 12.11 it is clear that nothing changes if 
functions ip E S are replaced by ipt G * (smoothness of 92 is not used 
there) . □ 

3.3 Proof of Propositions 12.51 and 12.61 

As Zt (and g are centered Gaussian, it suffices to prove convergence of 
the covariances. By (|2.5p and (|2.1U|) . we have for ip,ip £ S 

E{Z T ,(p)(Z T ,Tp) = (E9)I T + (Var 9)II T , (3.13) 

where 

It = \ I St<p{x)StiP(x)cIx — J St<p(x)T2StiP(x)(Ix 
F T \Jr Jr 

IIt = / STf(x)T2STil>(x)da 
Jr 

Observe that St<p{x) = T(p{Tx) (see (|2.9|) ) and recall that Tt(f(x) = e~ t],x],a (p{x) 
Hence, by Plancherel's identity, after an obvious substitution, we have 

T 1 



tx. 



I T = -3— {I- e- 2 W aT - a Mx)iP(x)dx. 
Ff 2tt Jr 

This and (l3~T3]) imply that if Var 9 = and F T = \/T l ~ a , then 

pin r 

lim E{Z T ,v)(Z T ,^) = — / <p(x)j>(x)\x\ a dx. (3.14) 

On the other hand, if Var 9 > and Ft = yT, then lim^oo It = and the 
covariance functional of Zt tends to Var 9-^ L tp(x)ip(x)dx = Var J R ^(^^(a^dx. 
This proves Proposition 12.51 

Furthermore, if Var 9 = 0, a < 1, and \& = {ipt)t>o have one of the forms 
(jl.6p - (jl.8p . then it is easy to see that the same argument can be repeated 
with (p = ipt, ip = ifis (and Ft = VT 1_a ) and we obtain 

1 E9 I " ~ 

7pi=z E QTtQTa = — / Mx)i) s {x)\x\ a dx (3.15) 

by (|3TT4D . 

It is not hard to calculate the right hand side of (|3.15|) explicitely. De- 
pending on the specific form of (tpt)t>o it is, up to a multiplicative constant, 
the same as the right hand side of (jl.ip . (|1.2|) . (|1.3p . respectively. To derive 
this, it is convenient to use the formula 

ii \x\ 

valid for a < 1. For example, for ipt of the form (|1.7|) (recall that ^ = 
![o,t] - J-H.o]) w e have 
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/ 



It xTTT, itv, , /" (1 - cos(xs))(l - cos(xt)) 
MxWs(x) \x \ a dx = 4 / : 7 _ — —dx 

Jr fI a 

(1 — cos(xs)) + (1 — cos(xt)) — i(l — cos(a;(s — i)) — i(l — cos(rc(s + t)) 
i m 



4a 



; ( S 1 - a + t 1 - Q -^(|s-i| 1 - Q + ( S + t) 



1-a- 



(cf. (fL2]) ). Note that this is a part of the proof of (pT9|) . □ 
3.4 Proof of Proposition 12.71 

Denote (p a {x) = ip(x — a). (a) follows immediately from (|2.5p since 71(</3 a ) = 
(7I(/?) a , and for $ of the form ([0]) ^ t+a - ^ s+a = ip? ~ V's • 

(b) is a consequence of (|2.7p , (|2.8|) and of the fact that under assumptions 
on * 

(tp v (x) - ifj u (x))T2{i>t - ip s ){x)dx > 0. 



This is obvious if t i— > ipt(%) is nondecreasing for each a; S R, and if this 
property holds on M + and ipt is odd for i > 0, then, by symmetry of p2 we 
can write this integral as 

/■OO /'OO 

2/ / {p2{x-y) - P2{x + y)){ipt{y) - 4> s (y))(4>v(x) - *fj u (x))dydx, 
Jo Jo 

which is non-negative since p2 is unimodal. 

(c) follows from (|2.5p . ()2.8|) (y>i</?2 = if <^i,^2 have compact supports 
and r is large) using the well known fact that \\m T ^ 00 T l+a p2(x) = c uni- 
formly in x on compact sets. 

Finally, to obtain (d), it suffices to observe that 



E(X(l),(p) z < C / <p 2 {x)dx, 
Jr 

since J R P2(x)dx < oo. □ 
3.5 Proof of Theorem I2T91 

For brevity, we will assume that at the beginning each interval [j,j + 1), j 6 
Z, contains exactly one particle. We then modify the notation of Section 
13. H writing instead of i?- 2 ' 1 , pi instead of p\\i an d h(x) = p^ — x instead 
of hi t i. 

Again, we use the central limit theorem, showing first the convergence 
of covariances and then applying the Lyapunov criterion. 

(a) Fix ip,ip £ S. Note that by the self-similarity property of the a-stable 

motion, has the same law as $ T _ a . Then, similarly as in 13.11 (c.f. ()3.3j) . 
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(|3.5p - (|3.8p ). we can write 

E(Z T ,ip) (Z T ,ip) 

— E E [<P& + + ^) - rjtt E + ^ + r ) 



r l- 



j62 



,x + h(x) 



1 



where 



Eip( — — + y)Eij){ + z)p T - a (y)p T - a (z)dydzdx 



(3.16) 



v?(|; + y)i>{^z + y)p T -« (y)dydx 



<f(j; + y)i>(j; + z)pr-« (y)Pr-°> (z)dydzdx 



(3.17) 



-E 



,x + h(x) . lf x + h(x) . ,x .,.x . 

v( — + i/M — 7^ + y)- <p( f + y)il>(f + y) 



p T - a (y)dydx, 



(3.18) 



jr- 1 + V)m ^ +z)- ip(- + y)^{- + z) 

p T - a {y)p T - a {z)dydzdx. (3.19) 



In Bt we substitute x' = % and pass to Fourier transforms. Then 



1 

2^ 



B T = T a — I [l- e -^ x \ a 



and it is clear that 



As if G 5, we have 



lim Bt = — 

T->oo 7T 



(p(x)ijj(x) \x\ a dx. 



,x + h(x) . ,x , 

<p( — 7^ ± + y)-^Y + y) 



c 
< — 

~ T 



by the Lipschitz property and (|3.4p . 

As before, denote <j){x) = . Using 

\ip(u + v)\ < C(l + u 2 )4>(v), u,v G 

we get 

l^(^ ± ^ + y)l<2^(| + y ) 

for T > 1, again by 



(3.20) 



(3.21) 



(3.22) 
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\U T \ < 



()3,2ip and ()3.22p imply, after the usual substitution, 

c r , , c 



T i- 



{(p(x+y)+\ip{x+y)\)p T - a {y)dydx 



T i- 



(<f>(x)+\if>(x)\)dx, 



\V T \ < 



C 



J^l—a 

c 



{<f>{x + y) + \ip(x + z)\)p T - a (y)p T - a (z)dydzdx 

3 

(4>(x) + \ip(x)\)dx. 



Hence 



lim Ut = lim Vt = 0, 

T-s-oo T^oo 



and by (|3.16p and ()3,20p we obtain 

lim E{Zt,(p){Zt,iP) 



ip(x)ip(x)\x\ a dx. 



(3.23) 



(3.24) 



Now, fix m such that (1/2)(1 — a)m > 1. Similarly as in the proof of 
Theorem O (see (EHUD . (T5TTTT) . (j5TT^|> ). to finish the proof it suffices to 
show that for any ip £ S, 



lim DrOp) = 0, 

T— >-oo 



(3.25) 



where 



1 ^ I J 



1 



(1 — a)m 

T 2 



i x + /i(x) 
£M Vy) p T - a {y)dydx 



(cf. ipnBjl ). ([53ZD implies that 

C 



D T (<p) < 



(l — a)m 



- [ (P m (x)dx 
1 Jr 



as T — > oo, 



by our assumption on m. 

(b) The argument of part (a) can be repeated. The only difference is 
that the functions of the forms (|1.6p - (|1.8p are not Lipschitz, so (|3.2ip does 
not hold for them. Nevertheless, instead of (|3.2ip . (|3.22p we can use the 
estimates 



x 



,x + h(x) . _ „ 

i[o, t] ( T w +y) < i^ijCy+v), 

as well as similar inqualities for lr_ tj0 ]- This permits to derive analogues of 
(fCTD . ([H^SD and d22S]). We omit details. □ 
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3.6 Proof of Propositions I27TT1 and I27T41 

Proposition 12.111 can be obtained by the same argument as that applied to 
derive Theorem 2.1(a) in [5], therefore we omit the proof. We will sketch 
the proof of Proposition 12.141 

Fix a > 1 and consider tpt = l[o,t] — If— t,ol- For arbitrary cti, . . . , a n € R, 
ti, ■ ■ ■ ,t n > denote 



fc=i 



It suffices to show that 



lim £ e W)W>> 
r-s>oo 



e 2 



(3.26) 



(3.27) 



where -fC is, up to a constant, the covariance function of nsfBm (see (|1.2p ). 
As if) is an odd function, by (jl.lOp and the Poisson initial condition we have 



cxp 



i-hp Jo i>(*+#r)dr 



Ee Vt 



1 dx 



(3.28) 



where $ is the standard a-stable process. 
We start with calculating the covariance 

1 

f 
2 
T 
1 



E I — I ip t (x + $ r )dr / ip s {x + -d u )dudx 
o Jo 

ipt (^Tu-rtps (x)dxdudr 

V>t (x) e~ n ' x ' " V> s (x) dxdudr 



o 



T /-T-r 

Jo 



4 /* /■ (1 - cos tx) (1 - cos «g) / _ ^ ^ 
1 r (1 — cos tx) (1 — cos sec) 



4 

T->oo 7T 



|2+o 



-dxdr = K(t, s) 



(3.29) 



(cf. (|2.19p ). Since ip is odd and •& is symmetric, (|3.27[) will follow from (|3.29p 
and ()3.28p if we show that for any t > 



lim R T = 0, 

T->oc 



(3.30) 



where 



Rt = T2 \ E 



T \ 4 

?/>i(x + # r )cZr ) g?x. 



o 



Similarly as before, 



I-RtI < 



< 



c_ 

2^2 



Ci 



T i-T r T r T 



Jri JT2 Jrz 

My ~ z ) 



-(r-2-ri)\x\ c 



M x - y) 



-(r8— ra) 



dxdydzdr^dr^dr2dr\ 



1 — cos tx 1 — cos t{x — y) 1 — cos i(y — z) 1 — cos tz 1 — e 



F - 2/1 



I?/ - z| 



M 1 2/ 1 



Fl 12/ 1 

dxdydz. 
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Fix any < e < ~. Using (1 — e ) \y\ a < T 1 £ \y\ a£ it is easy to see 
that 

. . C2 /* 1 — cos tx 1 — cos , 

Uvr < — — / i-: n dxaz, 

1 1 — Te / „ 1 il+a 1 1 1+a ' 

hence ([330]) follows. □ 
Proof of Theorem 12.151 

The idea of the proof is similar to the one used in the proof of Proposition 
12. 141 Again, we take ip of the form (|3.26[) and want to prove (|3.27p . Let 
N X)T denote the empirical process of the system started from one particle 
at point x. 

By the Poisson initial condition 

Ee i{Y T (im =exp |y (Ee^fi^'^-AllTdxX (3.31) 

We will need the following facts about N X,T 
E(N*> T ,p) =T r <p{x) (3.32) 

E (N?' T , Vl ) (N*' T , if 2 ) = T r {ipiT u -rV2)(x) 

pr 

+ V T T v ((Xr-v¥i){T u -vV2))(z)dv, r<u (3.33) 



E (N?' T , ip) 3 = T r (<p 3 ) + 3V T [ %-u {9 v ,t(u, -)T u <p) (x)dx, (3.34) 

Jo 



where 

g v ,T(u,x) = E(N?' T ,^) 2 

(see [20] and pi] (A.4.3) and (A.4.4)). 

(|3.32p and the fact that ip is odd imply that 



E I IT I ( N ?' T > ^) drH T dx = 0. 

Jr F T Jo 



(|3.27p will follow from (|3.3ip if we prove that 



A T (s,t):=^E f [ [ (N^ T ,i, s ){N^ T ^ t )drdudx — > K(s,t) 

■f T JM. Jo JO T^oo 

(3.35) 



T 

and 



B T (t):=^E[ [ (N?' T ,ip t )dr 
Jti Jo 



dx — ► (3.36) 

T— >oo 



By (|3.33p and the form of Ft, 



T r T 



+Vr J ^ s (x)(T u +r-2v^t){x)dv J dxdrdu 
=h{T) + I 2 (T), (3.37) 
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where 

f-T 



h{T) = 1 -^ F III ^ s {x)M^ u - rMa dxdrdu 

l-Kl V T Jo Jo Jr 

4 f (1 — costx)(l — cos sx) f 1 1 — e~ Tr l x l 



■drdx 



kV t Jr x 1 Jo \x\ 

< g ^'^ -)• 0, as T — ^ co, (3.38) 



and 



T /-T 



i" 2 (T) = — / $ 8 {x)Mx) / / / e-(" + "- 2f -)l^ Q dr^dx 

Ml J R Jo Jv Jv 

2 f (l-coste)(l-cossx) Z" 1 (1 - e - Tv N Q ) 2 , , 

-dvdx 



n Jr x 2 Jo \x\ 2a 

2 f (l-cos^)(l-cos.x) ^ (3 3g) 

7T ./to I.tI Q 



since a < |. 



(13351) now follows from ^37}-<^M and flffi). 

For brevity, denote / = \ipt\- By (|3.34p we have 



Br{t)<S^-[ ^ E(N^ T ,f) 3 drdx 
Jr Jo 



T 

--Ih(T) + II 2 {T) + Ih(T), (3.40) 



where 

1 rp 3 



Ih{T) = T 2 3 / / T r f(x)drdx = I f{x)dx -»■ (3.41) 

as T — > 00 (see (gJS)| ), 



T 



HxVt 




II 2 (T)=SJ— — / / / T u f\x)T u f{x)dudrdx 



JO 



-n/),/-^T^o, (3.42) 



(T u f(x)) / T v (Tu-vf) 2 (x)dvdudrdx 
Jo 




as T — >• 00, 

1 

7 I 

<d(/)^1^0 (3.43) 
#| 

as T ^ 00 by ipnB]L Thus (pM) is a consequence of M|l - (pll3l . This 
completes the proof. □ 
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